arXiv:hep-ph/0211234v2 24 Jan 2003 


Axial charge segregation and CP violation dnring the electroweak phase transition 

with hypermagnetic fields 

Alejandro Ayala'^’, Jaime Besprosvanyt 
]Instituto de Ciencias Nucleares, Universidad Nacional Autonoma de Mexico, 

Apartado Postal 70-543, Mexico Distrito Federal 04510, Mexico. 

^Instituto de Fisica, Universidad Nacional Autonoma de Mexico 
Apartado Postal 20-364, Mexico Distrito Federal 01000, Mexico. 

We study the scattering of fermions off the true vacuum bubbles nucleated during a first order 
electroweak phase transition in the presence of a background hypermagnetic field. By considering 
the limit of an infinitely thin bubble wall, we derive and analytically solve the Dirac equation in 
three spatial dimensions for left and right chirality modes and compute reflection and transmission 
coefficients for the case when the fermions move from the symmetric toward the broken symmetry 
phase. Given the chiral nature of the fermion coupling with the background field in the symmetric 
phase, an axial asymmetry is generated during the scattering processes. This mechanism is thus 
shown to provide an additional source of CP violation within the standard model. We discuss 
possible implications for baryon number generation. 
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I. INTRODUCTION 

It is well known that the standard model (SM) of elec¬ 
troweak interactions meets all the requirements -known 
as Sakharov conditions 0- to generate a baryon asym¬ 
metry during the electroweak phase transition (EWPT), 
provided that this last be of first order. As appealing 
as this scenario might be, it is also well known that nei¬ 
ther the amount of CP violation within the minimal SM 
nor the strength of the EWPT are enough to generate a 
sizable baryon number Hi 

Nonetheless, it has been recently shown that in the 
presence of large scale primordial magnetic fields it is 
possible to generate a stronger first order EWPT iilS 
The situation is similar to what happens to a type I 
superconductor where an external magnetic field mod¬ 
ifies the nature of the superconducting phase transition, 
changing it from second to first order due to the Meissner 
effect. 

The presence of magnetic fields during the evolution of 
the early universe has neither been proved nor has it been 
ruled out. Temperature anisotropies measured by COBE 
place an upper bound Bq ~ 10“® G for homogeneous 
fields {Bq refers to the intensity that the field would have 
today under the assumption of adiabatic decay due to 
the Hubble expansion) Q. In the case of inhomogeneous 
fields their effect must be searched for in the Doppler 
peaks and in the polarization of the cosmic microwave 
background radiation Q. The future satellite missions 
MAP and PLANCK may reach the required sensitivity 
for the detection of these last signals. 

Independently of their origin, there are several interest¬ 
ing cosmological consequences of the possible existence 
of primordial magnetic fields such as the generation of 
galactic seed fields and their influence on big bang nucle¬ 
osynthesis M- 

Recall that for temperatures above the EWPT, the 


SU(2)xU(I)v symmetry is restored and the propagat¬ 
ing, non-screened vector modes that represent a mag¬ 
netic field correspond to the U(l)y group instead of to 
the U(l)em group, and are therefore properly called hy¬ 
permagnetic fields. 

In previous works El , it has been shown that another 
interesting cosmological consequence of the presence of 
such fields during the EWPT is the production of an ax¬ 
ial charge segregation in the scattering of fermions off the 
true vacuum bubbles within the SM. This charge segre¬ 
gation is a consequence of the chiral nature of the fermion 
coupling to hypermagnetic fields in the symmetric phase. 
The analysis relied on the solution of an effectively one¬ 
dimensional Dirac equation for fermions moving perpen¬ 
dicular to the bubble wall. 

In this paper, we lift this restriction and solve the Dirac 
equation in three spatial dimensions for the two chirality 
fermion modes in the presence of an external hypermag¬ 
netic field. Working in the limit of an infinitely thin bub¬ 
ble wall, we find analytical solutions -valid for fermion 
wave functions described by small quantum numbers- 
from where we compute reflection and transmission coef¬ 
ficients and show explicitly that these differ for left and 
right-handed chirality modes incident from the symmet¬ 
ric phase. Since these coefficients are related to the ones 
obtained by studying fermion incidence from the broken 
symmetry phase by CPT and unitarity, we show how an 
axial charge segregation is built during the process. Such 
axial charge segregation provides a bias for baryon over 
antibaryon production and thus a seed for the non-local 
baryogenesis scenario MMM We also discuss how 
this mechanism is equivalent to an extra source of CP 
violation. 

The outline of this work is as follows: In Sec. im we 
write the Dirac equation for the left and right-handed 
chirality modes propagating in a background hypermag¬ 
netic field during the EWPT where the spatial profile 
of the Higgs field is taken in the limit of an infinitely 


2 


thin wall. In Sec. cni we find the three dimensional so¬ 
lutions to this equation and discuss their properties. In 
Sec. ITVI we use the solutions to compute reflection and 
transmission coefficients for axial fermion modes moving 
from the symmetric phase toward the broken symmetry 
phase. Finally in Sec. El we conclude by looking out 
at the possible implications of such axially asymmetric 
fermion reflection and transmission. 


II. DIRAC EQUATION FOR CHIRAL 
FERMIONS IN A BACKGROUND 
HYPERMAGNETIC FIELD 

During a first order EWPT, the properties of the nu¬ 
cleated bubbles depend on the effective, finite tempera¬ 
ture Higgs potential. In the thin wall approximation and 
under the assumption that the phase transition happens 
when the energy densities of both phases are degener¬ 
ate, it is possible to find a one-dimensional analytical 
solution for the Higgs field spatial profile (j){z) called the 
kink. This is given by 

^( 2 ;) ^ 1 -I- tanh(z/A), (1) 

where 2 : is the coordinate along the direction of the phase 
change and A is the width of the wall. 

When A is small compared to the mean free path, scat¬ 
tering within the wall is not affected by diffusion and the 
problem of fermion reflection and transmission through 
the wall can be cast in terms of solving the Dirac equation 
with a position dependent fermion mass, proportional to 
the Higgs field [ 13 , where the proportionality constant 
is given by the particle’s Yukawa coupling. In physical 
terms, the wall can be considered thin when the spatial 
region where the fermion mass changes is small compared 
to other relevant length scales such as the particle mean 
free path. 

In order to obtain analytical solutions, let us further 
simplify the problem considering the limit of an infinitely 
thin bubble wall, namely A —> 0. In this case, the kink 
solution becomes a step function, 8 ( 2 ;), and consequently, 
the expression for the particle’s mass becomes 

m{z) = moO(z). ( 2 ) 

For 2 : < 0 particles move in the region outside the bub¬ 
ble, that is, in the symmetric phase where particles are 


massless. Conversely, for 2 ; > 0, particles move inside the 
bubble, which represents the broken phase where parti¬ 
cles have acquired a finite mass mo- 

In the presence of an external magnetic field, we need 
to consider that fermion modes couple differently to the 
field in the broken and the symmetric phases. 

For 2 ; < 0, the coupling is chiral. Let 

= ^(1 + 75)^ 

= ^( 1 - 75 )^ (3) 

represent, as usual, the right and left-handed chirality 
modes for the spinor dt, respectively. Then, the equa¬ 
tions of motion for these modes, as derived from the elec- 
troweak interaction Lagrangian, are 

(i^ - - m(z)'i>B = 0 

(i^ - = 0, (4) 

where ya.L are the right and left-handed hypercharges 
corresponding to the given fermion, respectively, g' the 
U{1)y coupling constant and we take = (0, A) repre¬ 
senting a, not yet specified, four-vector potential having 
non-zero components only for its transverse spatial part, 
in the rest frame of the wall. This last requirement means 
that the corresponding magnetic field will point along the 
direction perpendicular to the wall. 

The gauge field profile should correspond to our choice 
of Higgs field profile and thus the change in the magnetic 
field strength should take place in the same small spa¬ 
tial region over which the change in the Higgs field takes 
place. An exact treatment that considered a continuous 
change in the Higgs field would require that the gauge 
field be also continuous across the interface. Nonetheless, 
since we are pursuing insight from analytical results, we 
consider the simplified situation in which the gauge field 
is constant on either side of the wall and discontinuous 
at the planar phase boundary. A discussion on the con¬ 
sequences of considering a continuo us g auge field profile 
can be found in the second of Refs. pT|. 

The set of Eqs. can be written as a single equation 
for the spinor 'I' = d>fl-|-4'j;,by adding them up 


-4 ^9'4+ 14 + ^9'4- - to ( 2 :)|«'= 0 . 


( 5 ) 


r 


Hereafter, we explicitly work in the chiral representation of the gamma matrices where 

0 cr 


7° = ( _? ^ ) 7 = 


-cr 0 


75 = 


I 0 
0 -I 


■ ( 6 ) 
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Within this representation, we can write Eq. as 


Since = 0, 


|z(^ - - to(z)|^' = 0 , 

where we have introduced the matrix 


g = 


^g'l 0 \ 

0 ^g'l ) ■ 


(7) 

( 8 ) 


We now look at the corresponding equation in the broken 
symmetry phase. For z > 0 the coupling of the fermion 
with the external field is through the electric charge e 
and thus, the equation of motion is simply the Dirac 
equation describing an electrically charged fermion in a 
background magnetic field, namely, 


= 0 . 


(9) 


In the following section, we explicitly construct the solu¬ 
tions to Eqs. 0 and m with a constant magnetic field, 
requiring that these match at the interface z = 0 . 


III. SOLVING THE DIRAC EQUATION 


= Pf = P, . 


(14) 


In this representation ' 01^2 and ^34 correspond to compo¬ 
nents with right- and left-handed chirality, respectively. 
Notice that, since in the symmetric phase fermions are 
massless, the differential Eqs. ca do not mix compo¬ 
nents with different chirality. 

Let us now work in cylindrical coordinates and look for 
solutions of the form 


4'(r, (l),z) 


/ \ 

/ 2 (r)e +*'^/2 

/3(r)e-*^/2 

V /4(r)e+*^/2 / 




(15) 


Using the explicit components of the vector A(r) that 
produce a constant magnetic field along the z axis, it is 
easy to show that 0 


For definiteness, let us consider a constant magnetic 
field B = Bz pointing along the z direction. We first 
look for the stationary solutions of Eq. 0 corresponding 
to fermions moving in the symmetric phase. Writing 

^(r,t) = e-*^‘^(r), ( 10 ) 


where 


4'(r) 


/ \ 
Mr) I 
Mr) 

\ Mr) ) 


( 11 ) 


and working in the chiral representation of the gamma 
matrices, Eq. 0 becomes explicitly 

E'ljjiir) - {P^ - iP^)'iP2{r) - PzMr) = 0 

EMr) - {Px + iPy)Mr) + PzMr) = 0 

EMr) + {Px - 'iPy)Mr) + PzMr) = 0 

EMr)+ {PMiPy)Mr)-PzMr) = 0 , ( 12 ) 

where we introduced the definitions 


± ( A ± i|. T 7^.1 , (16) 


which, together with the definitions 


= 

pR,L ^ 


yR,Lg, 


B 




n 

n. 


R,L _ 
1 — 

R;L _ 
2 — 


,L^l/2 ^ 


M-Lp^ 

(^«Up«U)i /2 ^2 


d 


dp 


,R,L 


- 1 - 


+ 1 + 


(l-l) 


rfR,L 


pR.L 


,(17) 


make the set of Eqs. 03 to be explicitly written as 


{E-k)^M^)+^P2f2{p^) = 0 
{E + k)h{p^)+i'R^h{p^) = 0 
{E+ k)h{p^) = 0 

{E-k)MP^)+^Pih{p^) = 0 . (18) 


'pR 

pi 


-iV 


-iV 


vM 


y^g' 


A. 


By acting on Eqs. itHHl with iPf, and 

respectively, we obtain the set of second order differential 
(13) equations 
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\ d{p^Y 2 4 4p^ J 

^ , Afl _ ([^ _ ^ _ Jll 

d{p^ Y 2 4 Ap^ j 

f,L d- , ^ P" a-l)M 

\ dp^ 2 4 4/?-^ / 

d(pi)2 ^ (ipi ^ 2 4 4pi J 


fiip^) 

f2{P^) 

h{p^) 

U{P^) 


0 

0 

0 

0 , 


(19) 


where 


A^'^ = 


E2 _ (fcfl,i)2 
^^R,L/ 


( 20 ) 


The eigenvalues corresponding to the requirement 
that the solutions vanish as p^'^ oo are such that 


/3(P^) = C'3/„-l,.(p^) 

/4(p^) = iCiIn,s{p^), (23) 

where C'i _....4 are constants and In,sip) are the Laguerre 
functions given in terms of the Laguerre polynomials 
Qr^ by 


X^'^=n. ( 21 ) 

n = / + s is called the principal quantum number which 
must be a non-negative integer n = 0,1,2,3... From 
Eq. 1201), the spectrum is given by 

E = ^(fc«T)2 (22) 

The solutions to Eqs. m are 

/l(p«) = Ciln-l,sip^) 
f2ip^) = iC2ln,sip^) 


satisfying the normalization condition 

nOO 

/ dp/2,(p) = l. (25) 

Jo 

Using the recurrence relations for the Laguerre polynomi¬ 
als together with the set of Eqs. 113, one readily obtains 
that the constants Ci ^2 and C' 3_4 are not independent and 
thus that the solution 4'(r) for fermions moving in the 
symmetric phase is given by 


J 


4'(r,^,z) = 




U3 




olK Z 


(26) 


We now turn to finding the solutions to Eq. ®, namely, 
for fermions moving in the broken symmetry phase. By 
a procedure similar to the above, it is easy to see that 
the solutions are given by ^ 


T(r, (t>,z) 


! Diln'-i,s'ip)e^P'-^^^ \ 

iD2ln',s'ip)e'^^'‘^ 

\lD4ln',s'ip)e'^^'‘^ J 


(27) 


where 


p = 


1 = ^B. (28) 

Notice that in the second of Eqs. we have used 

that since in the broken symmetry phase there should 
not be a propagating component corresponding to the 
field |j, Q , the parameter representing the magnetic 
field strength B' is related to B and Weinberg’s angle 9w 
by 


B' = 


B 

cos 9w ’ 


(29) 


which in turn implies that the coupling of the fermion 
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with the magnetic field in the broken symmetry phase 
eB' becomes 

eB' = g'B . (30) 

The spectrum is given hy E = + with 


n' = I' + s' being a non-negative integer n' = 0,1, 2 .... 
The constants are not independent from each 

other and are given in terms of the two constants Gi 
and G 2 by 


(D, 

D 2 

D 3 

V A 


/ 

= Gi 

V 




mo 

(E-k) 


mo 

0 

1 


In order to describe the motion of fermions in the whole 
domain, we should now proceed to join the solutions in 
Eqs. (EHll and (ITTIl across the interface at z = 0. To 
this end, we work in the approximation where the in¬ 
cident fermions from the symmetric phase make energy 
and angular momentum conserving transitions to states 
with the same principal quantum number in the broken 
symmetry phase, namely n = n' and 1 = 1', which re¬ 
quires s = s'. This approximation corresponds to the 
classical picture in which the effect of a sudden change 
in the coupling with the external field produces only a 
change in the radius of the orbit around the field lines. 
This is a reasonable approximation as long as the num¬ 
ber of quanta in the transverse modes for a given energy 
is not too large. This can be quantified by considering 

I 


\ 

/ 

(E+k) 

mo 

\ 

+ G 2 


yj — k'^ —rn^ 



rriQ 




1 


/ 

1 

0 

/ 


the overlap of the radial functions In,s and In',s' {p) 

for I = I' as a function of s' — s, as shown in Fig. 1. We 
find that our approximation is valid for low values of s 
and 1. Otherwise, interference with modes with addi¬ 
tional nodes associated with large values of these quan¬ 
tum numbers is not negligible, which means these modes 
need to be considered in the complete solution. 

For fermions incident from the symmetric phase, the 
complete wave function contains incoming as well as 
outgoing components, whereas in the broken symmetry 
phase the wave function corresponds to only outgoing 
waves. To describe the scattering of right- and left- 
handed chirality modes, let us consider them separately. 

For right-handed incident modes, continuity of the so¬ 
lution across the interface yields the system of conditions 


lE + k^ 
E-k^ 


-b 


'E-k^ , « 

E + kR ' 2 


. 

mo 


Gf- 


E + k" 
mo ^ 


Gi, 


l + G^ 


lE-k^ 
E + k^ 


G^ 


'E-k 
, mo 


nki _ 

Ly Q — ^9 


Gf. 


Gf- 


'iR I V E^ — k^ — iitQ j 

mo I ^ 


(32) 


On the other hand, for left-handed incident modes, conti¬ 
nuity of the solution across the interface yields the system 


of conditions 


J 


/ jE-k^ 

1 V E + k^ 


Gf 


f \/F2 -k^-mf, \ 

( — 


E + k 
mo 


G 


L 

2 
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( I E + k^ 
\\j E-k^ 


ck = - 


i + ck 


lE-k^ 
E + k^ 


= Gk 


Ck = G\ 


E-k 

mo 


Gk- 


yjE'^ - /t2 ^ 
mo 


Gk 


(33) 


Eqs. (T^ and (P|l completely determine the constants 
C' 2 ? 3 : and C'|; 3 , Gf 2 , respectively. 


IV. TRANSMISSION AND REFLECTION 
COEFFICIENTS 


The fact that the solutions to the systems of Eqs. 1211 ) 
and are not the same, means that there is the pos¬ 
sibility of building an axial asymmetry during the scat¬ 
tering of fermions off the wall. To quantify the asym¬ 
metry, we need to compute the corresponding reflection 
and transmission coefficients. These are built from the 
reflected, transmitted and incident currents of each type. 
Recall that for a given spinor wave function dt, the cur¬ 
rent normal to the wall is given by 


nOO n27r 

J= dp . (34) 

Jo Jo 

For right-handed incoming waves, the incident current 
is thus given by 


Jk 


2k^ 

E-k^ 


(35) 



FIG. 1: Overlap integral for the functions In,s{p^) and 
dn',B'{p) normalized to the maximum value in each case for 
/ = Z' as a function of s' — s for three different values of the 
pairs s and 1. Notice that when s and I are large, the overlap 
is maximum for s' s. 


whereas the reflected and transmitted currents 
are given respectively by 

2 fc« 


jR _ 
"'ref - 


E + k^ 


R\2 


(Ck) 


2k^ 


E + k^ 


{Ck) 


R\2 


jR _ 
"^tra 


y^E'^ — k'^ — m^ 


mo 


/ y'E^ - fc 2 _ , 
mo 

- (Gf)2 + (Gf)^ 


Gf 


Gk 


kE + k^ 

Gk 

\ rno J 


kE-k^ 

Gk 

\ rno J 



(36) 


On the other hand, for left-handed incoming waves, the 
incident current is given by 


7- 

U-iy- 


2k^ 

E-k^ 


(37) 


and the reflected and transmitted currents are 

given respectively by 

^ {ck)^ ^ 


jL _ 

"'ref - 


E + k^^ 


jL _ 
"'tra ~ 


- fc2 _ . 

mo 

' sjE"^ - : 


\ mo 

- {Gk)^ + {Gkf. 


E + k^ 


Gk 


L\2 


{Ck) 


Gk 


kE + k^ 

Ck 

\ rno J 


I'E-k's 

1- 

\ rno J 



(38) 


Recall that in the symmetric phase, where fermions are 
massless, the chirality and helicity operators commute 
and thus chirality modes are also eigenfunctions of helic¬ 
ity. From Eq. and the representation of the gamma 
matrices, Eq. we see that for right-(left) handed chi¬ 
rality modes, the large components correspond to right- 
(left) handed helicity modes. Since scattering of the 
wall does not change the spin direction of the imping¬ 
ing particle, right-(left) handed helicity modes reflect as 
left-(right) handed and transmit as right-(left) handed 
modes. To emphasize this point, for the reflection and 
transmission coefficients, let us denote by lower case let¬ 
ters r and I the right- and left-handed helicity modes. 
These coefficients are given as the ratios of the corre¬ 
sponding reflected and transmitted currents, to the inci¬ 
dent one, explicitly 

Rr^i = dkef/'^'mc 










































FIG. 2: Coefficients Ri^r and Rr^i as a function of the mag¬ 
netic field B scaled by rriQ for E/mo = 3, n = 1, and ya = 4/3, 
Vl ~ 1/3. The value for the U{1)y coupling constant is taken 
as g' = 0.344, corresponding to the EWPT epoch. 

= Jtm/4c- (39) 

The corresponding coefficients for the axially conjugate 
process are 

Rl^r = 

Ti^i = 4a/4c- (40) 

Using Eqs. (I35II to (13011 . together with the solutions to 
Eqs. (EU and E3, it is straightforward to show that 

Rr —4“ Rr — >r — 1 

Rl^r+Tl^l = 1 (41) 

Figure 2 shows the coefficients Ri^r and Rr^i as a func¬ 
tion of the magnetic field B scaled by TOq for E/mo = 3, 
n = 1 and hypercharge values for top quarks ya = 4/3, 
yL = 1/3, with g' = 0.344, as appropriate for the EWPT 
epoch. Notice that when 6^0, these coefficients ap¬ 
proach each other and that the difference grows with in¬ 
creasing field strength. 

Figure 3 shows the reflection and transmission coeffi¬ 
cients as a function of energy scaled by the fermion mass 
toq. Figure 3a shows the coefficients Rr^i and Tr^r and 
Fig. 3b the coefficients Ri^r and Ti^i for B/m^ = 0.5 
and n = 1, j/fl = 4/3, j/l = 1/3, g’ = 0.344. Since the 
solutions in Eqs. ini) are computed assuming that the 
transmitted waves are not exponentially damped, their 
energies have to be taken such that E > yJm^E dny. 
For these energy values, there is thus no need to consider 
the contribution from negative-energy solutions Ea 

V. DISCUSSION AND CONCLUSIONS 

In this paper we have derived and solved the Dirac 
equation describing fermions scattering off a first order 
EWPT bubble wall in three spatial dimensions, in the 
approximation of an infinitely thin wall, in the presence 
of a magnetic field directed along the fermion direction 
of motion. In the symmetric phase, the fermions couple 




FIG. 3: Reflection and transmission coefficients as a function 
of energy scaled by the fermion mass for B/m/o = 0.5, n = 1, 
yn = 4/3, yL = 1/3, g' = 0.344. (a) Shows the coefficients 
for incident, right-handed helicity modes and (b) for incident, 
left-handed helicity modes. 


chirally to the magnetic field, which receives the name of 
hypermagnetic, given that it belongs to the U(1)y group. 
We have shown that the chiral nature of this coupling 
implies that it is possible to build an axial asymmetry 
during the scattering of fermions off the wall. We have 
computed reflection and transmission coefficients show¬ 
ing explicitly that they differ for left and right-handed 
incident particles from the symmetric phase. The results 
of this calculation are in agreement with those previously 
found in Refs. 0 , where the motion of the fermions was 
effectively treated as one-dimensional. 

The three dimensional treatment of the problem allows 
for a clearer physical picture of the fermion motion in the 
presence of the external field: Suppose that the original 
fermion motion is not parallel to the direction of the field 
lines and therefore that its velocity vector contains a com¬ 
ponent perpendicular to the field. In this case, due to the 
Lorentz force, the particle circles around the field lines 
maintaining its velocity along the direction of the field. 
The motion of the particle is thus described as an over¬ 
all displacement along the field lines superimposed to a 
circular motion around these lines. These circles corre¬ 
spond to the different quantum levels, labeled by n, the 
principal quantum number. We see that the originally 









different angles of incidence all result in the same overall 
direction of incidence. Also, since the relevant quanti¬ 
ties that determine the strength of the coupling are the 
magnitudes of the parameters 7 '^’'^ and 7 , the generated 
axial asymmetry is independent of whether the fermion 
moves parallel or anti-parallel to the field lines. 

What is the relation of this axial asymmetry to CP 
violation? Recall that in the SM, CP is violated in the 
quark sector through the mixing between different weak 
interaction eigenstates to form states with definite mass. 
However, in the present scenario, taking place within the 
SM, no such mixing occurs since we are concerned only 
with the evolution of a single quark (for instance, the top 
quark) species. The relation is thus to be found in the dy¬ 
namics of the scattering process itself and becomes clear 
once we notice that this can be thought of as describing 
the mixing of two levels, right- and left-handed quarks 
coupled to an external hypermagnetic field. When the 
two chirality modes interact only with the external field, 
they evolve separately, as described by Eq. ll^ . It is 
only the scattering with the bubble wall what allows a 
finite transition probability for one mode to become the 
other. Since the modes are coupled differently to the ex¬ 
ternal field, these probabilities are different and give rise 
to the axial asymmetry. CP is violated in the process 
because, though C is conserved, P is violated and thus 
is CP. 

We also emphasize that, under the very general as¬ 
sumptions of CPT invariance and unitarity, the total ax¬ 
ial asymmetry (which includes contributions both from 
particles and antiparticles) is quantified in terms of the 
particle (axial) asymmetry. Let pi represent the number 
density for species i. The net densities in left-handed 
and right-handed axial charges are obtained by taking 
the differences Pl — Pi and pa — pp, respectively. It is 
straightforward to show H that CPT invariance and 
unitarity imply that the above net densities are given by 

PL-Pl = if - - Rl-,r) 

PR- PR = ir - f){Rl^r - Rr^l) , (42) 

where /® and are the statistical distributions for par¬ 
ticles or antiparticles (since the chemical potentials are 
assumed to be zero or small compared to the temper¬ 
ature, these distributions are the same for particles or 


antiparticles) in the symmetric and the broken symme¬ 
try phases, respectively. From Eq. (03), the asymmetry 
in the axial charge density is finally given by 

{pL - pl) - iPR - Pr) = - f){Rr^l - Rl^r).m 

This asymmetry, built on either side of the wall, is disso¬ 
ciated from non-conserving baryon number processes and 
can subsequently be converted to baryon number in the 
broken symmetry phase where sphaleron induced transi¬ 
tions are taking place with a large rate. This mechanism 
receives the name of non-local baryogenesis [I1II1II1II3 
and, in the absence of the external field, it can only be re¬ 
alized in extensions of the SM where a source of CP viola¬ 
tion is introduced ad hoc into a complex, space-dependent 
phase of the Higgs field during the development of the 
EWPT [ 3 . 

An interesting possibility to extend the scope of the 
present work is to study the scattering of fermions in 
the presence of topologically non-trivial configurations of 
hypermagnetic fields such as the so called, hypermagnetic 
knots which themselves can seed the ba^'on asymmetry 
of the universe in extensions of the SM . 

Since another consequence of the existence of an ex¬ 
ternal magnetic field is the lowering of the barrier be¬ 
tween topologically inequivalent vacua [^, due to the 
sphaleron dipole moment, the use of the mechanism dis¬ 
cussed in this work to possibly generate a baryon asym¬ 
metry is not as straightforward. Nonetheless, if such pri¬ 
mordial fields indeed existed during the EWPT epoch 
and the phase transition was first order, as is the case, 
for instance, in minimal extensions of the SM, the mech¬ 
anism advocated in this work has to be considered as 
acting in the same manner as a source oi CP violation 
that can have important consequences for the generation 
of a baryon number. 
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